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ABSTRACT 

It is known that theory of MOND with spherical symmetry cannot account for the 
convergence K-map of Bullet Cluster 1E0657-558. In this paper, we try to set up a 
Finslerian MOND, a generalization of MOND in Finsler spacetime. We use Ric = 
to obtain the gravitational vacuum field equation in a four-dimensional Finsler space- 
time. To leading order in the post-Newtonian approximation, we obtain the explicit 
form of the Finslerian line element. It is simply the Schwarzschild's metric except for 
the Finslerian rescaling coefficient f(v) of the radial coordinate r, i.e. R = f(v)r. 
By setting f(v) = yl — (GMao/w 4 ), we obtain the famous MOND in a Finslerian 
framework. Taking a dipole and a quadrupole term into consideration, we give the 
convergence n in gravitational lensing astrophysics in our model. Numerical analysis 
shows that our prediction is to a certain extent in agreement with the observations of 
Bullet Cluster 1E0657-558. With the theoretical temperature T taking the observed 
value 14.8 keV, the mass density profile of the main cluster obtained in our model is 
the same order as that given by the best-fit King /3-model. 
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1 INTRODUCTION 



It has long been noticed that according to Newton 's inverse-square law of gr avity, the observed baryonic matter cannot provide 
enough force to attract the matter of the galaxies ( OortjB"932l ; Zwickv.lll933h . This inconsistency has been confirmed by a large 
numb er of observation in the past thirty years, to name a few, the ve l ocity dispersions of dwarf Spheroidal galaxies (|VogtJ 
1995) and the flat rotation curves of spiral galaxies II Rubin et al .. 1980; IWalter et al.Jl20oih . et al.. Postulatin g that galaxies 



are s urrounded by massive, non-luminous dark matter is the most widely adopted way to solve the problem l|de Blok et al 



2008). The dark matter hypothesis has dominated astronomy and cosmology for almost 80 years. However, up to now, no 
direct observations have been firmly tested. 

Some models have been built as an alternative of the dark matter hypothesis. Their main ideas are to assume that 
the Newtonian gravity or Newton's dynamics is invalid on galactic scales. The most successful and famous model is MOND 
l Mikrom.|[l983l ). It assumes that the Newtonian dynamics does not hold on galactic scales. The MOND paradigm is based 
on the following assumptions: (i) It introduces a new physical constant ao = 1-2 x 10~ 8 cm/s 2 . (ii) The law of gravity returns 
to Newton's gravity while ao —¥ 0. (iii) The law of gravity is given as om = yGMao/r in the deep-MOND limit, ao — > oo. 
As a phenomenological model, MOND explains well the flat rotation curves of thousands of spiral galaxies with a simple 
formula and a univers al constant. In particu lar, it naturally gives the well-known global scaling relation for spiral galaxies, the 
Tully-Fisher relation (jTullv fc Fisher.lll977h . The Tully-Fisher relation is an empirical relation between the total luminosity 



of a galaxy and the maximum rotational speed. It is of the form L oc wj^ax, where a m 4, if the luminosity is measured in the 
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2000) showed that the Tully-Fisher relation appears to be convergent 



near-infrared region. Tully and Pierce ( Tullv fc PierceJ 
in the near-infrared region. McGaugh ( Mc GaughJl2005r ) investigated the Tully-Fisher relation for a large sample of galaxies, 
and concluded that the Tully- Fisher relation is a fundamental relation between the total baryonic mass and the rotational 
speed. MOND (|MilgromJll98a ) predicted that the rotational speed of galaxy has an asymptotic value v 4 ^^^ = GMao, which 

explains the Tully-Fisher relation. 

By introducing several scalar, vector and tensor fields, Bekenstein (IBekenstein .1120041 ) rewrote the MOND into a covari- 
ant formulism (TeVeS). He showed that the MOND satisfies all four classical tests of Einstein's general relativity in Solar 
system. How ever, MOND still faces challenges. T he strong and weak gravitational lensing observations of Bullet Cluster 
1E0657-558 ddowe. Randall fc MarkevitchJ 2007 ) cannot be explained by MOND and its Bekenstein's relativistic version 



I Angus et al.. 12006. 2007 ). The ICM (intracluster medium) gas accounts for most of the Bullet Cluster's mass. Clowe et al. 



Clowe. Randall fc MarkevitchJ 120071 ) had r econstructed the surface ma ss density Ti(x,y) from the Chandra space satellite 



X-ray image of the ICM gas. Moffat et al. (jBrownstein fc MoffatJ 120071 ) had shown that the E-map of the ICM gas of the 
main cluster can be well fitted with a King /3-model density profile. The King /3-mode l is a radial distribution of the mass 
density for a nearly isothermal and isotropic gas sphere. On the other hand, Clowe et al. ( Clowe. Randall fc MarkevitchJl2007V ) 
had reconstructed the convergence re-map from the strong and weak gravitational lensing survey. The re-map indicates that 
additional gravitational force is needed for explaining the Bullet Cluster. The center of gravitational force deviates from the 
center of the ICM gas. And the distribution of gravitational force does not possess spherical symmetry. Most of theories of 
modified gravity, such as Bekenstein's relativistic version of MON D, only consider radial force. Also, most of the mass density 
profile of dark matter, such as the NFW profile 1 Navarro] 2007 ). only contrive radial (isotropic) distributions. All of them 
cannot explain the observations of the Bullet Cluster. 

The distribution of gravitational force in Bullet Cluster is anisotropic. To describe anisotropic force, one should introduce 
the multipole fields. The dipole contribution vanishes if one takes the center of ICM gas as the coordinate origin. Monopole 
contribution plus quadrupole contribution are needed to account for the observations of Bullet Cluster. In fact, Milgrom 
gave a Quasi-linear formulation of MOND (QUMOND)l Milgrom. 2010T ). which involves the quadrupole contribution. Usually, 
MOND effects vanish in Newtonian regime. However, Milgrom showed that quadrupole effect appears even in hig h-acceleration 



syste ms. Besides, Angus et al. presented an A r -body code for solving the modified Poisson equation of QUMOND (lAngus et al. 



2012). They used it to compute rotation curves for a sample of five spiral galaxies from the THINGS sample ijWalter et al. 



2008) and concluded that taking gas scale-heights of the gas-rich dwarf spiral galaxies (and stellar scale-heights of stellar 



dominated galaxies) as free parameters is vital to make precise conclusions about MOND. Other interesting results were also 
obtained in their work. 

On the other hand, besides Bekenstei n's TeVeS, there are other 'MONDian' theories (for example, the Einstein-aether 
theory i Zlosnik. Ferreira fc Starkmanl2o"o7 )). Both the Bekenstein's and the Einstein-aether theory admit a preferred reference 
frame and broken local Lorentz invariance. It can be reasonably inferred that the local Lorentz invariance violation (LIV) is 
an intrinsic feature of MOND. If this is acknowledged, there follows a conclusion: the space structure near a galaxy is not 
Minkowskian even at long distances from the galaxy center. It depends on the rotational velocity of the galaxy considering 
the relationship between the Tully-Fisher relation and MOND. 

Finsler gravity, which is based on Finsler geometry, has the feature s mentioned above. Thus it is natural to postulate that 
Finsler gravity is a covariant formulism of MOND. Finsler geometry ( Bao. Chern fc ShenJ 2000h is a natural generalization 
of Riemannian geometry with the latter as a special case. The length element of an arc on a Finslerian manifold depends not 
only on position but also on velocity, which induces anisotropy. Preservation of the fundamental principles is a prerequisite for 
Finsler gravity as well as the results of gen eral relativity. A new geometry (i.e. the Finsler geometry) involves new spacetime 
symmetries. Kostelecky ( KosteleckvJ 2011 ) has shown that LIV is closely related to the Finslerian geometry. Effective field 
theories have been studied in his paper for explicit LIV effects in Finslerian spacetime. 

In l|Li fc Changj|201ll ). we presented the vacuum field equa tion in Finsler gravity, and have given the solution of field 
equation under weak field approximation. In ( Li fc Chang J 2012 ). we presented the Newtonian limit in Finsler gravity and a 
covariant formulism of MOND. We studied the spacetime structure of MOND with properti es of Tully-Fisher r elation and 
Lorentz invariance violation. A Finsler spacetime has less symmetry than a Minkowski one ( Li fc ChangJ I2OI0I ) . Multipole 
effects such as dipole and quadrupole contributions, which embody space anisotropies, should be considered in Finsler gravity. 
In this paper, we try to construct a Finslerian MOND, a generalization of MOND in Finsler spacetime, and use it to explain 
the observations of Bullet Cluster. 

The rest of paper are organized as follows. Section 2 is dedicated to the theory used for the numerical analysis and is 
separated into five parts: Section 2.1 is about the basic concepts of Finsler geometry; in Section 2.2, we give a brief introduction 
of the special relativity in Finsler spactime; in section 2.3, we extend Pirani's argument to a Finsler spacetime to obtain the 
gravitational vacuum field equation in Finsler gravity; in Section 2.4, the Newtonian limit in Finsler spacetime is presented; 
in Section 2.5, under post-Newtonian approximation we give the Finsler structure. In Section 3, we give the convergence re 
in our model. Section 4 is about the numerical analysis which contains two parts: in Section 4.1, we consider the dipole and 
quadrupole contributions to the Finslerian MOND in the calculation of the convergence re of the Bullet Cluster; in Section 
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4.2, we obtain the mass density of the main cluster given the observed value of its surface temperature and compare it to the 
best-fit King /3-model. Numerical results are presented in 2D as well as in 3D figures. Conclusions and necessary discussions 
are presented in Section 5. Demonstrations of certain approximations in Section 3 is presented in the Appendix. 



2 FORMULISM OF FINSLER GRAVITY 
2.1 Basic Concepts 

Finsler geometry is based on the so called Finsler structure F. F is a non-negative real function which has the property 
F(x, Xy) = \F(x,y) for all A > 0, where x represents position and y = 4^ represents velocity. The fundamental tensor is 
given as (|Bao. Chern fc Shen]|200cl) 

d d 



9»v 



r V2 J 



dy^ dy v 

The arc length in Finsler space is given as 

dr 1 dr n 
plx 1 ... x n - _ . . . zr_W r 



(1) 



(2) 



A more detailed discussion about F can be found in Section 5. Hereafter, we adopt the following index gymnastics: Greek 

indices in lower case run from 1 to 4, while Latin indices in lower case (except the alp habet n) run from 1 to 3. 

The parallel t ransport has been studied in the framework of Cartan connection ( MatsumotoJ Il98fj| ; lAntonelli fc Rutz 



198rj; li zabo, 2008). The notation of parallel transport on a Finsler manifold mea ns that the length F (^) is constant. The 
geodesic equation on a Finslerian manifold is given as (|Bao. Chern fc Shen,ll200oh 

+ 2G M = 0, (3) 



dr 2 
where 



G M = -a 
4 a 



d 2 F 2 A 

-,y 



of 2 

dx" 



are called the geodesic spray coefficients, r is arc length on the Finsler manifold. Obviously, if F is a Riemannian metric, then 

(5) 

where *y** \ is the Riemannian ChristofTel symbol. Since the geodesic equation @ is directly derived from the integral length 



G = ip v \v v > 



F 



dx \ 
d-r) 



dr, 



(6) 



the inner product ( ^/ g M „ = F (^) j of two parallel transported vectors is preserved. 

On a Finslerian manifold, there exists a linear connection - the Chern connection I Chern J H) IS). It is of torsion freeness 



and almost metric-compatible 

-pa a o 

1 T p.V Q 



N P N p N p 

A " A * _L A P 

Ax^p — —j^- + A vX p —p- 



(7) 



where 7 a MJ/ is the formal ChristofTel symbols of the second kind with the same form of Riemannian connection. iVt is defined 
as ATt = 7ta£/ a — A ll vX 'f X a py a 'y /3 and Axp,u = -j -^jtt{F 2 ) is the Cartan tensor (regarded as a measurement of deviation 

from the Riemannian Manifold). In terms of the Chern connection, the curvature of Finsler space is given as 



5r A 



8xv 



8x" 



1 a/j 1 k,v 



y X K.ll 



(8) 



where A = t^tt - N"-^. 

oxP oxt 1 H- ay' 



2.2 The Null Set F = and Finslerian Special Relativity 

In Finsler geometry, the Finsler structure F is defined as a non- negative C°° function on the entire slit tangent bundle TM\0, 
i.e. F : TM — > [0, oo). It ensures that the integral length ((2| always makes sense (since a negative arc length is not acceptable 
in mathematics). In physics, for a gravity theory, the quantity F 2 dr 2 represents the line element of spacetime (which is also 
called 'proper time interval' in some references). A positive, zero and negative F correspond to time-like, light-like ('null') 
and space-like curves respectively. For massless particles, the stipulation is F — 0. 
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One should notice that many Finslerian geometric objects like Ricci scalar involves the Finsler structure F. It might 
be invalid to describe the massless particles at first glance. However, the ambiguities caused by F — can be removed 
by re-parameterizing the formulae with some other parameter a such that F(cr) 7^ 0. The property of Finsler structure 
F(x, Xy) = XF(x, y) guarantees that the length L is independent of the choice of curve parameter. Under a given parameter 
change r = C(a), £ > 0, the length L is of the form L(r) = £ F (x, |f ^) dr = £ F (x, |f) da = L(a), where r and a 
are both curve parameters and y = dx/dr (or y = d x/da). The same trick has been played in general relativity for massless 
particles which has g M „ d f d f = (|Weinbergjil"972l l. 



Einstein's special relativity contains two basic principles. The first one states that the law of nature are invariant under 
the Lorentz transformation. The other one states that the speed of light is the same for all inertial observers. In mathematical 
language, the first principle manifests that the length of four-dimensional Minkowski spacetime are invariant under the 
coordinate transformation. The way of describing spacetime symmetry in a covariant language (the symmetry should not 
depend on any particular choice of coordinate system) involves the concept of isometric transformations. In fact, the symmetry 
of spacetime is described by the so called 'isometric group'. The generators of isometric group is directly connected with the 
Killing vectors. The isometric group of Minkowski spac etime is the Loren tz group. In Finsler spacetime, the generation 
Minkowski spacetime is the locally Minkowski spacetime ( Li fc Chang. 20ld ). In order to investigate the isometric group of 



Finsler spacetime, we need to find its Killing vectors. 

It is convenient to discuss the isometric transformations under an infinitesimal coordinate transformation for x 

x"=x'" + eVi", (9) 

together with a corresponding transformation for y 

r = y» + ^y", (io) 

where |e| <C 1. Under the coordinate transformation ([9]) and (|10[) . to first order in ej, we obtain the expansion of the Finsler 
structure, 

dF 9U M dF 

my)=Fl X ,y) + eV»— + ey»— — , (11) 

where F(x,y) should equal to F(x,y). Under the transformation Q and (|10|) . a Finsler structure is called isometry if and 
only if 

F(x,y) = F(x,y). (12) 
Deducing from the (|11|) . we obtain the Killing equation Kv(F) in Finsler space 
dF OF 

K v (F) = V^+y-^— =0. (13) 
v ' dx» dx v dy^ 

Substituting F = ^/g^( x ,y)y' J 'y 1 ' into (|13[1 . the Killing equation Kv(F) could be written in terms of Finsler metric g til ,(x,y) 

cto. +9a »dxv +9 ^ l dx> 

Here, we must point out that the Killing equation (1141) is different from the one that obtained by Yano ( Yanol 1957 ). This 



+ + g Pl >%- ) y a y" = 0. (14) 



due to fact that the term ^rn-y v ^§^§-y a y^ in (|14[) equals 0. In general, a non-Riemannian Finsler spacetime admits less Killing 
vectors than a Minkowski spacetime. It is well known that a n-dimensional Minkowski spacetime has n(n+ l)/2 independent 
Killing vector s. In a n-dimensional Finsler spacetime, the maximal number of independent Killing vectors is n(n — l)/2 + 1 
( Szabo. 20081 ) . Our physical interest concentrates on four-dimensional spacetime. It means that a four-dimensional Finsler 



spacetime has seven independent Killing vectors at most. It is well k nown that the isom etric group is a Lie group on the 
Riemannian manifold. This fact also holds on the Finslerian manifold ( Deng fc HouJ 20021 ). 



The Finsler structure of locally Minkowski spacetime has the form F = W g^v{y)y IJ 'y l/ ■ The metric of locally Minkowski 
space is a function only of the fiber coordinate y. It could be regarded as a generation of Euclidean space or Minkowski 
spacetime. For it is a flat space in Finsler geometry. In a locally Minkowski spacetime, we find that 

V" = a M (15) 



is solutions of Killing equation (|13|) . where a M are constants. It means that a locally Minkowski spacetime has translation 
symmetry. Then, a four-dimensional locally Minkowski spacetime at most have a three dimensional subgroup of Lorentz group. 
The Lie algebra of locally Minkowski spacetime is also a subalgebra of Poincare algebra. Therefore, we complete the first step 
of constructing the counterpart of special relativity. 

The second principle of special relativity states that the speed of light is the same for all inertial observers. It also involves 
the concept of causal structure. In locally Minkowski spacetime, general speaking, the speed of light is not the same for all 
inertial observers. The coefficients of speed of light are constant but not equal to each other. Such property reflects that the 
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locally Minkowski spacetime is an anisotropic spacetime. For example, a locally Minkowski spacetime of Randers form is given 
as 

Fdr = yj dt 2 — dx\ — dx\ — dx\ — ^rdxi. (16) 

The null set of locally Minkowski spacetime is given by the equation Fdr — which describes massless particles like photons. 
It is obvious that the speed of light c — {1,1, \/f }■ 



2.3 Extension of Pirani's Arguments 

In this paper, we introduce the vacuum field equation in a way first discussed by Pirani (jPirani ] |l964lRutzJl998h . In Newton's 
theory of gravity, the equation of motion of a test particle is given as 



dr 2 



d.i 



(17) 



where (j> = </>( x ) i s the gravitational potential and rf 3 — diag(+l,+l,+l) is the Euclidean metric. For an infinitesimal trans- 
formation x % — ^ x l -f- e£ l (|e| 1), the equation (|17|l becomes, to first order of e, 



dr 2 dr 2 dx 1 dx Wx k 

Combining equations ()17[) and (|18|) . we obtain 

dr 2 V ? dx*dx k ~ ? k ' 

For the vacuum field equation, one has H\ — \/ 2 (j> = 0. 

In general relativity, the geodesic deviation gives a similar equation 

Dt 2 
where 



CR%, 



(18) 



(19) 



(20) 



„„^t—^t—- Here, R^„„ is the Riemannian curvature tensor. 'D' denotes the covariant derivative along the 

A up dr dr ' A up o 



curve a; M (f). The vacuum field equation in general relativity gives R^xu = ( Weinberg] 1972T I. This implies that the tensor 



is also traceless, R = R^ = 0. 

In Finsler spacetime, the geodesic deviation yields ( Bao. Chern fe Shen. 2000T I 



2tM 



D £ 
Dt 2 



(21) 



where R^ 



RJ 1 ^—^t-. Here, i?, M is Finsler curvature tensor defined in d8l), 'D' here denotes covariant derivative 

A up dr dr 'A up Y — V J 



■^p- = ^h- + r^, A (;r, ^p). Since the vacuum field equations of Newton's gravity and general relativity are of similar 

forms, we may assume that vacuum field equation in Finsler spacetime has similar requirements as in the case of Netwon's 
gravity and general relativity. It implies that the tensor in Finsler geodesic deviation equation should be traceless, = 0. 
In fac t, we have proved t hat the analogy of the geodesic deviation equation is valid at least in a Finsler spacetime of Berwald 

type ( Li fc Chang. 20081 ) . We assume that this analogy still holds its validity in a general Finsler spacetime. 

In Finsler geometry, there is a geometrical invariant — the Ricci scalar Ric . It is of the form ( Bao. Chern fc ShenJ 200oT ) 

dG" x <9 2 G" „„x <9 2 G M dG" <9G A> 



Ric = R^, 



1 

F 2 



' dx^ 



I) 



dx x dy^ 



+ 2G" 



dy x dy^ dy x dyi 1 



(22) 



The Ricci scalar depends only on the Finsler structure F and is insensitive to the connection. For a tangent plane FIcT^M 
and a non-zero vector y £ T X M, the flag curvature is defined as 

gxnR^u x 



K(n,y) = 



(23) 



F 2 g pg uPu e - (g aK yu K ) 2 ' 

where u £ II. The flag curvature is a geometrical invariant and a generalization of the sectional curvature in Riemannian 
geometry. The Ricci scalar Ric is the trace of i? M „, which is the predecessor of flag curvature. Thus the value of Ricci scalar 
Ric is invariant under the coordinate transformation. 

Furt hermore, the significance of the Ricci scalar Ric is very clear. It plays an important role in the geodesic deviation 
equation {hi fc ChangJl201ll , l2012llBao, Chern fc Shenj[200ol ). The vanishing of the Ricci scalar Ric implies that the geodesic 
rays are parallel to each other. It means that it is vacuum outside the gravitational source. 

Therefore, we have enough r easons to believe that the gravitational vacuum field equation in Finsler geometry has its 
essence in Ric = 0. Pfeifer el al. (|Pfeifer fc Wohlfarthj|2012l ) have constructed gravitational dynamics for Finsler spacetime 
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in terms of an action integral on the unit tangent bundle. The stipulation Ric — here is compatible with their results of 
gravitational field equatio: 



-gr; 
1 



2.4 The Newtonian Limit in Finsler Spacetime 

It is well known that the Minkowski spacetime is a trivial solution of the Einstein's vacuum field equation. In Finsler spacetime, 
the trivial solution of the vacuum field equation is called 'locally Minkowski spacetime'. A Finsler spacetime is called a locally 
Minkowski spacetime if there is a local coordinate system {a;' 1 }, with induced tangent space coordinates {y 11 }, such that F 
depends not on x but only on y. The locally Minkowski spacetime is a flat spacetime in Finsler geometry. Using the formula 
(|22p , on e can see that a locally M inkowski spacetime is a solution of Finslerian vacuum field equation. 



In (|Li fc Chang,ll201ll . I2012T ). we assumed that the metric is close to the a locally Minkowski one r)^ v {y), 



(24) 



considering that the gravitational field /i M „ is stationary (thus all time derivatives of fr M „ vanishes) and the particle is moving 
very slowly (i.e. GM/r <C 1). The lowering and raising of indices are carried out by ?7 M „ and its matrix inverse rf" v . We found 
from Ric 

d 2 h a f 



n 



dxidxi* " 0- (25) 

In general relativity, one uses post-Newtonian approximation to study the motion of particle ( Weinberg] 1972 ). Before 
studying the motion of particle in Finsler gravity, we must deal with the concept of energy-momentum tensor in Finsler space- 
time. It is well known that the energy-momentum tensor is conserved (in the sense of covariant differentiation) and symmetric 
in general relativity. How ever, this is not the case in Finsler gra vity. The energy-momentum tensor is symmetric if the angular 
momentum is cons erved JPubrovin. Fomenko fe Novikovll 19991 ). Generally, the symmetry of angular momentum is broken in 
Finsler spacetime ( Li fc Chang.l | 201C ). Thus, the energy- momentum tensor is not symmetric in Finsler gravity. Similar situa- 
tions appear in torsion gravity ^Hammond 2002) in Riemann-Cartan geometry. Also, to satisfy the con servation law, besides 



the R icci scalar, additional terms that represent the "torsion effect" are needed in the field equation (Pf eifer fc Wohlfarth 



2012t ). Although these "torsion" terms would cause a difficulty to understand Finsler gravity, they could fortunately be omit 
ted. The reason is that these "torsion" terms do not contribute to the geodesic deviation equation, which determines the 
motion of particles in Finsler geometry. Furthermore, we concentrate only on t he motion of partic l e with zero spin in a weak 
gravitational field. Therefore, with similar steps to deduce the equation ([25} in ( Li fc Chang]|201ll . 2012 ). and by making use 



of the post-Newtonian approximation, we obtain the gravitational field equation in Finsler gravity 



n 



(r a ,a - TjVafjT^ 



(26) 



hoo,h nn are terms of the same order as GM/r, and the corresponding component of the energy- momentum tensor is Too. 
Finsler gravity should reduce to general relativity, if the Finsler metric <? M „ reduces to a Riemannian one. Thus, we find from 
(1261 that 



ij d 2 hpo 
dx % dxi 
B 2 h 

%2 ''nn 

dx % dxi 



—8Tv F Gprioo, 
8tt f Gprj n „, 



(27) 
(28) 



where p = Tqo/t?oq is the energy density of the gravitational source. In Finsler spacetime, the space volume of r)^ v (y) 
i Bao. Chern fc Shenjlioool ) is different from the one in Euclidean space. We used hf in (I27I28[) to represent the difference, 
where 



7VF 



^fgdx 1 A dx 2 A dx 3 



(29) 



g = det{r)ij) is the determinant of r\ij. 'A' denotes the 'wedge product[f|. The solution of (|27I28[) is given as 



The gravitational vacuum field equation given in llPfeifer fc Wohlfarthj2012l ) is g F ab d a d b K- -^ I K+2g F ab (V a S b +S a S b +d a VS 6 ) 



: 0. 



The Sa-terms can be written as S a = ^P^ba' wnere ^ d = an d ^dba are * ne coefficients of the cross basis dxA^p- iBao. Chern fc Shen.1 
2000). Considering that TL = R a ab y b = -R a dab V d y b = F 2 (£ d Rf ab ^ b ) = F 2 (g ab R ab ) = F 2 Ric and dropping the S a -terms (see the 
discussions about the 'torsion' terms in next section), one can see that Ric = is one of the solutions of the above equation. 

2 The derivations in the rest of this and the next subsections, if not specifically pointed out, are accurate to the first-order of hfi V . 

3 In Subsection 2.2, we have discussed the local symmetry of Finsler spacetime. It manifests that the symmetry of locally Minkowski 
spacetime is different from the Minkowski spacetime. The space length that determined by the symmetry of locally Minkowski spacetime 
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2GM 2GM 

noo = =— »7oo, = — =—r} nn , (SO) 

ix H 

where R 2 = r/ijx'x-' . In Newton's limit, the geodesic equation (|3]) reduces to 

<fx° n 0i dh 00 dx° dx° 
~dr 2 T dx i ~~dT~dT 
dV _ rf J dh Q0 dx° dx° 
dr 2 ~ ~2~~dxT~dT~dV 



= 0, (31) 
= 0. (32) 



The equation (|31[) implies that ^5- is a function of ftoo- Since |/iooj <C 1, ^f— could be treated as a constant in equation ([32 
Then, we find from (J32) that 

d 2 x i GM x* 



dx° 2 ~ R 2 R ' (33) 

where dx° 2 = rioodx dx°. The formula (J33J implies that the law of gravity in Finsler spacetime is similar to that in Newton's 
case. The difference is that the spatial distance is now Finslerian. It is what we expect from Finslerian gravity, because the 
length difference is one of the major attributes of Finsler geometry as compared to the Riemannian geometry. 

2.5 Finslerian MOND 

In l|Li fc Chang]|2012h . we have shown that Finsler gravity reduces to MOND, if the spacial part of the locally Minkowski 



metric of galaxies is of the form 

where ao = 1.2 x 10 _10 m/s 2 is the constant of MOND. In Finsler spacetime, the speed of particle is given as v' = ^ s = 
The radial coordinate in the locally Minkowski space-time of galaxies (|34[) can be written as 



yfq^xi = rJl - {^^) = rf(v), 



R= v/%iW=rWl- j-^ =rf{v), (35) 



where r — 8ijX % x 3 and v = 5ijV*v 3 . Substituting (|35|) back into (|33|l . we obtain the result of MOND 
r \ rao J 



where (J,(x) = x/\/x 2 + 1 is the interpolating function in MOND. 

In this paper, we try to consider multipole effects of Finslerian MOND, and use them to explain the observed K-map of 
Bullet Cluster. The Finslerian radial coordinate has the form R = rf(v). And without losing any generality in our discussion 
of the motion of particle in Finsler spacetime, we set r/oo to 1. Then, we obtain the Finsler structure in the post-Newtonian 
approximation from (1301) 



F W = (l-™M) dr 2 -(l+ 2 -^) dR 2 ■ ^ 2 



R J \ R J 

To first order in h, the geodesic spray coefficients are 



ir(dr + sin Qdtf). (37) 



G" = 



Iff y -—y /j. (38) 

Then, given ((38]), one can solve the geodesic equation ((3|. And by making use of the stipulation F = in ()37[) for photons, 
one could obtain the formula of gravitational deflection of light in Finsler spacetime. We skip the conventional calculations 
here. In fact, the Finslerian line element (|3TI) is simply the Schwarzschild's one except for a rescaling of the Euclidean radial 
coordinate r. It is also true for the geodesic equation. Thus, the deflection angle in Finsler gravity is a rescaling of Einstein's 
one. It is of the form 

a F = — — , R m = r m }(v m ), (39) 

lira 

where v m is the fiber coordinate for corresponding r m . Rm is the closest distance of the light path to the gravitational source, 
where r m is that in a Euclidean space. 

is also different from the Euclidean length. So does the unit circle a nd its related quantity-7r. H ere, we denote the Finslerian 7r by np. 
'A' is the 'wedge product'. For more details please refer to the book l lChern, Chen fc Lamj|2006h . 
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3 CONVERGENCE k IN FINSLER GRAVITY 

The ICM gas of Bullet Cluster can be well described by the King /3-model ( Cavaliere fc Femiano, 19761 ). Its mass density 



distribution is given as 

2\ "3/3/2 



P(r) = PO 1 + -o 



The surface mass density is given by integrating p(r) of equation (|40|l along the line of sight 



p(x 1 ,X2,x 3 )dx 3 



(40) 



(41) 



In th e limit 2:3"* 3> r c , the surface mass density is of the form (detailed discussions can be found in (jBrownstein fc Moffat 
2007|)) 

/ f 2\ "(3/3-l)/2 

E(0 = E 1 + K 



(42) 

where £ 2 = x\ + x\ is defined in the lens plane and Eo = y/Tvpor c V( 3|fl 2 ~ 1 ) /T(^). Moff at et al. For the observed IC M gas 
profile of the main cluster, the best- fit parameters of the King /3-model (|40f) are given as i Brownstein fc Moffat J 2007 ) 



P = 0.803 ± 0.013, r c = 278.0 ± 6.8 kpc, Po = 3.34 x 10 5 Mo/kpc 3 



(43) 



1976) 



In observations, the convergence kg measures the ratio of observed surface density to the critical surface density ijPeacock 



, , , , D L D LS _ E(xi,aj 3 ) 
K G = 47TG / — - p(xi,X2,x 3 )dx s = , 



D s 



(44) 



where E c = 



, j , , , , Ds is the distance between the source galaxy and the observer, Dl is the distance between the lens 
(Bullet Cluster) and the observer, and Dls is the distance betw een the source galaxy and the lens. For the Bullet Cluster, the 



critic al surface density E c takes a value of 3.1 x 10 9 Mo/kpc 2 ( Clowe. Gonzalez fc MarkevitchJ 2004 ; Brownstein fc Moffat 
2007h 



In general relativity, the convergence k maps the gravitational lensing effect of a given gravitational source. It can be 
expressed in terms of the deflection angle as 
1 D LS D L ^ 



D s 



(45) 



where £ = \/ x \ + x \- The deflection angle in Finsler spacetime was given as (|39[) . Substituting cyf into (|45p . we obtained 



1 D LS D L 
(tj? = — 

2 D s 



where 



kg 



1 D LS D L 



1 « r, 1 

-— V 5 a G +«gV 5 -— 
/(«) /(«) 



1 ID LS D L 1 

77T K G + =: Q GV 5 -— 

/(v) 2 D s f(v) 



(46) 



(47) 



~ KG ■ 



(48) 



2 

The first term is simply a rescaling of k given by general relativity. The second term depends on the specific form of /. It 
does not retain the linearity and the superposition principle of the point mass potential on mass m. It can also be neglected 
in the weak field approximation. We will demonstrate the second point in APPENDIX and show that for the two cases of / 
(see the next section) that were investigated in this paper, the second term is a few orders smaller than the first term and 
thus can be neglected. The Kp can be approximately given by 
1 

KF ~ -J7-, 
/(«) 

Here, we summarize the logic steps to deduce the convergence kf in Finsler gravity. First, we extended Pirani's argument 
of equation of motio n to the case of Finsler geo metry to get Ric = 0, which can be derived from an action integral on the 
unit tangent bundle (jPfeifer fc Wohlfarth.|[2012i ). Second, in post-Newtonian approximation we obtained gravitational field 
equation in Finsler gravity (|26[) . Third, in the Newtonian limit to first order of GM/R, we obtained the Finsler line element 
(|37|l . It is simply the Schwarzschild's metric except for the rescaling coefficient f(v) of the Euclidean radial coordinate. Then, 
we obtained the deflection angle (|39[) in Finsler gravity. Fourth, given the relation (|45p between the convergence k and the 
deflection angle a, we obtained the Finslerian convergence kf as (148(1 . Up till now, our formulae in Finsler gravity haven been 
presented on the tangent bundle. However, the physics of the astronomical observations lie in four-dimensional spacetime. 
We need a projection that translates the formulae on the tangent bundle into the ones on the manifold. Such a projection 
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stems from the solution of geodesic equation. The geodesic equation (|3} gives the relation y = ^ — y(x). It implies that f(v) 
could be written as a function of x by the relation y(x). Finally, after doing these steps, we obtain the Finslerian convergence 

k f ~ —3— k g , (49) 

where g(x) = f(y(x)). Given the surface mass density profile (|42p , we could obtain the numerical results of convergence K-map 
from the equation (l49ll . 



4 NUMERICAL ANALYSIS 
4.1 The Convergence K-Map 

The surface mass density profile (|42|1 with best-fit parameters (|43|l are shown in Figure [T] The main X-ray cluster is set at 
£ = kpc and the subcluster (the peak of which lies at £ ~ 400 kpc) is neglected in doing the best-fit with the King /3-model. 
The surface mass density (I42|l for the main cluster includes most of the ICM gas. It implies that the ICM gas profile of the 
Bullet Cluster is in approximate spherical symmetry. We will use the surface mass density (|42[) to calculate the convergence 

In this paper, our motivation is to construct a MOND-like theory in Finsler gravity, and use it to explain the observations 
of the Bullet Cluster. As mentioned in the introduction, a modified gravity theory is taken as a theory of MOND so long as it 
reduces to Newton's gravity while the MONDian constant ao — > and the Tully-Fisher relation holds for deep-MOND limit, 
aa — > oo. 

First, we propose a Finslerian MOND with spherical sy mmetry. The geode sic equation gives an approximate relation 



between the velocity and the modified gravitational potential 1 Li &: ChangJ 2012 ) 
3M _ GM 
~R~ ~~ rf(v 



2 G M GM 

v = S- = 7773- ( 5 °) 



+ VGMao . (52) 



If 

9M(r) = f{v{r))= + I . (51) 

we find from (|50[) that 

2 _ GM 
r 

It is a MOND theory with spherical symmetry. It should be noticed that the three-dimensional radial distance r equals the 
two-dimensional radial distance £ if one deals with the physics in the lens plane. So, in this section, we use r to represent the 
radial distance on the lens plane. Substituting (|51|l into (|49|l . we obtain the convergence K-map given by MOND theory with 
spherical symmetry. The result is shown in Figure [2] In Figure [2] one can find that a MOND theory with spherical symmetry 
cannot account for the reconstructed convergence K-map of Bullet Cluster. The convergence K-map of Bullet Cluster shows 
that the distribution of gravitational force is anisotropic. To describe the anisotropic force, we should introduce multipole 
fields. The dipole contribution vanishes if one t akes the center o f ICM gas as the coordinate origin. In fact, Milgrom gives a 
quasi- linear formulation of MOND (QUMOND) (jMilgrom Jlipioh , which involves the quadrupole contribution. 



Here comes our second step. We take the quadrupole effect into consideration in Finslerian MOND in a way that the 
quadrupole contribution appears even at large scales. The Finslerian parameter go:(r,9) now takes the form 



-1, m , . aor 2 ( GMao 2 Q , , s 
9q {r,$) = l + d [} + —frT- cos <? ex P(-r/c) 



(53) 



where the parameters b — 458 km/s and c = 220 kpc. In order to keep the Tully-Fisher relation, an exponential term exp(— r/c) 
is needed in (|53[) . Substituting (|53|) into (|49[) . we obtain the convergence K-map giving by MOND theory with monopole 
contribution plus quadrupole contribution. The result is shown in Figure [3] The monopole contribution plus quadrupole 
contribution can account for the main feature of the convergence K-map of Bullet Cluster, except for the asymmetry between 
the convergence of the main cluster and the subcluster. Until now, we only consider the effect of the spherical part (i.e. the 
main cluster of the S-map) of ICM gas. The dipole contribution vanishes as we take the coordinate origin to be the center of 
ICM gas. 

Here comes our final step, we consider the subcluster of the E-map and regard it as a perturbation. Equivalently, it could 
be regarded as a dipole contribution. Then, The Finslerian parameter g(r, 9) is of the form 



1 ,,, , a r 2 ( V GMa GMao 2 n , , \ ,r A \ 
9Q D { r ,0) = 1 + \j I H ^2 cos6 l exp(— r/c) H cos 0exp(-r/c) j , (54) 
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Figu re 1. The scaled S-map from X-ray imaging observations of the Bullet Cluster 1E0657-558, November 15, 2006 data release 
llClowe. Randall fc Markevitch.ll2006l . I2007T) . The peak of the main cluster is taken to be the referential center of the system, i.e. r = . 
The peak of the subcluster locates at £ ~ 398 kpc. A cross-section of the observed E-map, on a straight-line connecting the peak of the 
main cluster to that of the subcluster, is shown in solid black dots. The best-fit King /3-model for the surface mass density is shown in 
solid red. Negative radii have no particular meanings but only denote the left-hand-side of the map relative to the origin r = 0. They 
have the same physical interpretation as the positive- half radii after being placed an absolute sign '| |'. 



where the parameter a = 2b ~ 916 km/s. In formula (154} . the dipole term v a2 — ~ 1 for r ~ 780 kpc. It means that the 
dipole term in (|54[) becomes dominant at r ~ 780 kpc. Distance at this far almost reaches the boundary of the Bullet Cluster 
system. And it is suppressed by the exponential term exp(— r/c). Thus, it is justified to regard the dipole term in (|54[) as a 
perturbation. Substituting (|54[) into (|49[) . we obtain the convergence k which is predicted by the Finslerian MOND theory 
with the contribution of a monopole, a quadrupole and that of a dipole perturbation. The result is also shown in Figure [3] One 
can see the asymmetry between the convergence k peak of the main cluster and the subcluster, with the center of convergence 
kf for the system lying at a few kpcs away from the origin due to the dipole effect. The 3D figure is shown in Figure [4] and 
the observational data of the convergence k of Bullet Cluster is presented in Figure [5] for comparison. 

In the last section, we concluded that for our discussions, the equation (|46[) can be well approximated by equation (|48l) . 
Given the ([51} and 1(53)1. we check this point carefully in APPENDIX. 



4.2 The Isothermal Spherical Mass Profile 

At last, we give a discussion about the isothermal temperature of the ICM gas profile of the main cluster in Finsler gravity. 
The ICM gas profile of the main cluster can be regarded as a spherical and isotropic system. If it is in hydrostatic equilibrium, 
it satisfies the collisionless Boltzmann equation 

a(r)EE _^ = I^, (55) 
dr p dr 

where $ is the gravitational potential and ay is the velocity dispersion. Assuming an isothermal gas profile, ay is related to 
the isothermal temperature T of the ICM gas as 

i kT 



fim p 



(56) 



where k is the Boltzmann constant, fi ~ 0.609 is the mean atomic weight and m p is the proton mass. Substituting both 
formula (|56[1 and the density distribution (|40|l of King /3-model into ((55} , we obtain that 
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Figure 2. The K-ma p reconstructed from the strong and weak g ravitational lensing survey of the Bullet Cluster 1E0657-558, November 
15, 2006 data release l[cl owe. Randall fc Markevitch ] |200dl2007h . The solid black dots denote the cross-section of the scaled S-map from 
the X-ray imaging observations that presented in Figure [T] A section of the reconstructed re-map, on a straight-line connecting the peak 
of the main cluster to that of the subcluster, is shown in black stars. The peak of the main cluster locates at £ ~ —180 kpc and that of 
the subcluster locates at £ ~ 522 kpc. The £ = point is chosen to be the same with that of the S-map in Figure [T] The convergence Kp 
predicted by the Finslerian MOND is shown in solid green, for which the Finslerian factor jj^(r) is of the form {5TJ. The interpretation 
of negative radii is the same as that in Figure [l] 
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Figure 3. A cross-section on a straight line connecting the two peaks of the reconstructed re-ma p reconstructed from the strong and 
weak gravitational lensing survey of Bullet Cluster 1E0657-558, November 15, 2006 data release ijClowe. Randall fc Markevitch^. 2006, 
2007). The solid black dots denote the cross-section of the scaled E-map from the X-ray imaging observations that presented in Figure 
[T] Quadrupole contribution is considered in the Finslerian MOND in the prediction of convergence Kp, which is shown in solid red. The 
result that taking the subcluster's contribution of ICM gas into account as well as the quadrupole contribution is shown in solid blue, for 
which the Finslerian factor jQfl(r) is of the form Fill) . The interpretation of negative radii is the same as that in Figure [T] in addition 
that now the negative-half correspond to 9 = while the positive-half require that 8 = it. 
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Figure 5. The 3D K- map reconstructed from the strong and weak gravitational lensing survey of the Bullet Cluster 1E0657-558, November 
15, 2006 data release (Clowe. Ra ndall fc Markevitch] 2006. 2007). The interpretation of negative radii is the same as that in Figure [l] 
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Figure 6. The mass profile given by the best-fit King /3-model of the maincluster is shown in solid red. The mass profile derived from 
Finslerian MOND with quadrupole effect is shown in solid black. The results are presented on a logarithmic scale for both the M- and 
r-axis. 



Markevitch et al. (Markevitc h et al .. 2002) have presented the experimental value of the isothermal temperature of the 
mam cluster T = U.8±l 7 keV with 4 - 5% 

error. By making use of (|57|l . one can find that the Newton's gravitational force 
ajv = —GM/r cannot provide enough force to maintain the hydrostatic equilibrium. In Finsler gravity, the gravitational 
acceleration law is of the form 
GM _ GM 
' R 2 



Bf = 



(58) 



(rg(r)) 2 ' 

We neglect the dipole perturbation in our study of the isothermal temperature of the main cluster and only consider the 
quadrupole contribution in Finslerian MOND. Even this, the gravitational system is no more isotropic. Nevertheless, we could 
take the average of the radial force by integrating g(r,8) of ([53} over 8 



g( r Y 



1 

2^ 



r sMr .^ 1+ ^( 1+ ^« p( _, /c) ) 



(59) 



and use it to qualitatively study the hydrostatic equilibrium of an isotropic system. Substituting the g(r) of equation ([59 
into (15811. we obtain that 



GMa 




exp(— r/c) j 



(60) 



Her e, we take the temperat ure in formula (|57[) to be 14.8 keV, which is the experimental mean value given by Markevitch et 
al. ( Markevitch et al.. 20021 ). Then, by identifying the equation l|57[) with (I60[) . we obtain the mass profile of the main cluster 
of ICM gas in Finslerian MOND. In Figure [6] we compare it with the result of the best-fit King /3-model. It is shown that the 
two mass profiles have th e same order. It means that the Finslerian MOND with quadrupole effect agree with the observations 
[ Markevitch et al.Jliool ). 



5 CONCLUSIONS AND DISCUSSIONS 

In this paper, we try to setup a Finslerian MOND, a generalization of MOND in Finsler spacetime. We extended Pirani's 
argument to get the stipulation Ric = 0, from which we obtained the gravitational vacuum field equation in Finsler spacetime. 
Considering the correspondence with the post-Newtonian limit of general relativity, we got the explicit form of the Finslerian 
line element. It was simply the Schwarzschild's metric except for the Finslerian rescaling coefficient f(v) of the radial coordinate 
r. Given that f(v) = (GMao/v 4 ), we recovered the famous MOND in a Finslerian framework. By introducing a 
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quadrupole and a dipole perturbation term into the Finslerian MOND, we calculated the convergence k in gravitational 
lensing astrophysics. A qualitative-level numerical analysis showed that our prediction is in agreement with the observed 
rc-map of Bullet Cluster 1E0657-558. Given the observed value 14.8 keV of the isothermal temperature of the main cluster in 
our model, the predicted mass density profile of the main cluster is the same order as that given by the best-fit King /3-model. 

However, one should notice that the factor f(v) (i.e. g(r)) is determined by the local spacetime symmetry, which cannot 
be deduced from the gravity theory. It is not the fruit but a prior stipulation of the theory. The logic is: given a specific /(«), 
we then proceed to calculate the convergence K-map and the temperature of the main cluster. The coefficient f(v) in our 
model comes directly from the flat Finsler spacetime r]fj, v (y) (I24|l . In fact, while the Euclidean radial distance r — > oo, the 
Finslerian length element (|37[) reduces to 



dR -R\d6 z + sin Odtf). (61) 

It is simply the line element of flat Finsler spacetime ^^(y). The coefficient f(v) could be arbitrary in principle, since we 
suppose the metric is close to the flat Finsler spacetime rj^ (y) . Most of the galaxies could be regarded as a spherical system 
and described by a central modified gravitational potential. Therefore, to describe it in Finsler gravity, the Finsler parameter 
f(v) should be spherical. It means that the flat Finsler spacetime ^^(y) is not universal in cosmology. At present, the specific 
form of flat Finsler spacetime ^^(y) could be regard as an axiom in our theory of Finsler gravity. There is no physi cal 



equation or principle to constrain the form of it. Professor Shen(one of the authors of the book ( Bao. Chern fc ShenJ 20001 ) Vs 
description of Finsler geometry (private conversation) may help us in understanding what is a flat Finsler spacetime r),j, v (y) 
— "Riemann geometry is 'a white egg', for the tangent manifold at each point on the Riemannian manifold is isometric to a 
Minkowski spacetime. However, Finsler geometry is 'a colorful egg', for the tangent manifolds at different points of the Finsler 
manifold are not isometric to each other in general." In physics, it implies that our nature does not always prefer an isotropic 
gravitational force. It is also "colorful", as we have seen in case of Bullet Cluster 1E0657-558. 



APPENDIX 

We will show that for the two cases of / (or g) that were investigated in this paper, the second term in Eq.(46) are both a 

few orders smaller than the first term and thus can be neglected. 

-l 



Given that g M (r) = f(v(r)) = ( 1 + yj j for MOND, we get k f = j^kg + I^^V^ao- To get this 

result, we have replaced r with £ = \J x\ + x\. For a rough estimate of the magnitude of the second term, we write cea 
as olq = 4 Jj^ , where c ~ 3 x 10 s m/s = 9.71 x 10~ 12 kpc/s is the speed of light in vacuum. For the Bullet Cluster 
system, we have a total mass of M ~ 1O 14 M0 and a distance range of ^ £ ^ 1000 kpc. The constant for MOND is 
ao = 1.2 x 10~ 8 cm/s 2 ~ 3.84 x 10~ 30 kpc/s 2 . A simple arithmetic exercise shows that the magnitude of the second term in 
the expression of kf (i.e. the term | Dl ^^ l \J ~§M a G) is at about 10 -5 , which can be neglected comparing to the first term 
f(v(r)) KG ' °f wmcn value ranges from 0.1 ~ 0.4. 

The calculation in the last paragraph was carried assuming that M is a constant. If mass M is taken as a function of f, i.e. 



M = M(C), k f takes a form k f = j^kg + + { -%*"'■ ) «G^- ^/M. Using a G = ^f, the third 

-DlsDl \ V GMa . dM 



term on the right hand side of the above identity can be reduced into ( D d S s Dl ) a ° ■ ^fif /M. The term 4M. /M has an order 

dM 



of 10~ 3 kp^ 1 and y/GMa ~ 10~ 27 (kpc/s) 2 , giving that the whole third term ( ~ D ^ S ° L ) ^ G f a ° ■ yj§-/M ~ 10" 6 , which is 



negligible comparing to the first two terms in the expression of kf- Therefore, for <7m(^) = f(v(r)) = 1 + \ ^rW ) > ^ ne 



-i 



l 



convergence kf is given as kf — /(u(r)) KGm 

This conclusion still holds true for the anisotropic Finslerian MOND model we presented. Given that gQD{r,9)~ 



apr 2 ( 1 

gm y 1 



1 + cos exp(— r/c) + GM l a a cos 2 #exp(— r/c) J , we obtain the corresponding k as kqf, — g ^ fl ) K G + 



I Dl d s Dl a G^ 'g ^ 9 ) ■ We plot gq D (^,d) and V; gQ J'^ jpj as functions of £ respectively in Figure [7] One can see that 
Vg g() J~(£ g) — 10~ 2 kpc -1 . Again we can check that for the given parameters a = 2b ~ 916 km/s and c = 220 kpc, the term 

\ Dl d s Dl a G^£ gcjod 0) ~ 10~ 5 i s a ^ so small comparing to the first term jj^. e ^ kg- Thus it is justified to be neglected in 
the calculation of kf- 

Therefore, considering the above discussions, the convergence kf in our Finsler gravity model is given as kf — jnjhrsyj K G- 
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Figure 7. Plot of SqqC?,'?) and V vs. distance § in kpc. The blue solid line represents that for 6 = and the red one for 
8 = 7r. One can see that ^(9q d (S,S) — 10~ 2 kpc -1 . 
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